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ABSTRACT 



We have identified and quantified semi- major axis drifts in Near-Earth As- 
teroids (NEAs) by performing orbital fits to optical and radar astrometry of all 
numbered NEAs. We focus on a subset of 54 NEAs that exhibit some of the 
most reliable and strongest drift rates. Our selection criteria include a Yarkovsky 
sensitivity metric that quantifies the detectability of semi-major axis drift in any 
given data set, a signal-to-noise metric, and orbital coverage requirements. In 
42 cases, the observed drifts (~ 10" 3 AU/Myr) agree well with numerical esti- 
mates of Yarkovsky drifts. This agreement suggests that the Yarkovsky effect is 
the dominant non-gravitational process affecting these orbits, and allows us to 
derive constraints on asteroid physical properties. In 12 cases, the drifts exceed 
nominal Yarkovsky predictions, which could be due to inaccuracies in our knowl- 
edge of physical properties, faulty astrometry, or modeling errors. If these high 
rates cannot be ruled out by further observations or improvements in modeling, 
they would be indicative of the presence of an additional non-gravitational force, 
such as that resulting from a loss of mass of order a kilogram per second. We 
define the Yarkovsky efficiency fy as the ratio of the change in orbital energy 
to incident solar radiation energy, and we find that typical Yarkovsky efficiencies 
are ~10~ 5 . 

Subject headings: astrometry — minor planets, asteroids — minor planets, asteroids: 
individual (1999 RQ36, Aten, Apollo, Ganymed, Geographos, Hathor, Icarus, 
Orpheus, Ra-Shalom) — radiation mechanisms: thermal 
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Introduction 



Understanding how the Yarkovsky force modifies asteroid orbits has illuminated how 
asteroids and meteorites are transported to near-Earth space from the main belt and has 



allowed for deeper understanding of the structure of asteroid families ( IBottke et al. 



20061 ) 



The Yarkovsky force is necessary for acc urately predicting asteroid trajectories, including 



those of potentially haza rdous asteroids ( iGiorgini et al. 



2002 



2008 



Milani et all 120091 ). 



Chesley 



2006 



Giorgini et al. 



The Yarkovsky effect (or force) describes the process by which an asteroid's surface 
ther mal lag and rotation result in net thermal emission that is not aligned towards the 



Sun (IBottke et al 



2002b 



20061 ). The so-called diurnal component of the Yarkovsky effect 



operates as follows. A prograde-spinning object generally has a component of this surface 
thermal emission anti- aligned with the motion along the orbit, producing a net increase in 
the object's semi-major axis (i.e., da/dt > 0, where a is the semi-major axis). Conversely, a 
retrograde-spinning object generally has a component aligned with its velocity, shortening 
its semi-major axis (i.e., da/dt < 0). 

The maximum possible drift rate for any radiation-powered force acting on near-Earth 
asteroids (NEAs) can be obtained by equating the incident solar radiation energy in a given 
time interval to the change in orbital energy during the same interval. We find 

1 L Q 1 



da 3 

dt = * Y ~kH vT^e* GM Q Dp' 



1 



Goldreich fc Sari 



(120091 ). e is the 



where fy is an efficiency factor analogous to that used by 
eccentricity, L and M are the luminosity and mass of the Sun, G is the gravitational 
constant, and D and p are the effective diameter and bulk density of the asteroid. This 
equation exhibits the expected dependence on the asteroid area-to-mass ratio. In convenient 
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units, it reads 



da 



1.457 / fy \ fl km\ /1000 kg m 



10~ 3 AU/Myr. 



(2) 



dt v^^2 yiQ-sj \ D j v P 
Maximum efficiency (f Y =l) would convert all incoming solar radiation into a change in 
orbital energy. We will show in Section [3] that typical Yarkovsky efficiencies are fy ~ 10 -5 
and that typical rates are ~ 10~ 3 AU/Myr for kilometer-sized asteroids. The low efficiency 
and rates are due to the fact that it is the momentum of departing thermal photons that 
moves the asteroid. 



Chesley et al.l (120031 ) used precise radar ranging measurements to (6489) Golevka and 



reported the first detection of asteroidal Yarkovsky dri: 
da/dt = (-6.39 ± 0.44) x 10~ 4 AU/Myr Jchesley et al. 



t. The drift rate for this NEA of 



20081) corresponds to an efficiency 



f Y = 5x 1Q- 6 for £>=530 m and p = 2700 kg m" 



Vokrouhlicky et al. 



(120081 ) employed the Yarkovsky effect to link a 1950 observation 
to asteroid (152563) 1992 BF with a da/dt rate of (-10.7 ± 0.7) x 10" 4 AU/Myr. This 
corresponds to an efficiency fy = 7x 10~ 6 for D=420 m and p=2500 kg m -3 . If 1992 BF 
has a density closer to 1500 kg m~ 3 , the efficiency would be fy = 4 x 10~ 6 . 

There have been ot h er sea rches for the effects of non-gravitational forces in 



asteroid orbits. ISitarskil ( 119921 ) considered a semi-major axis drift in the orbit of 
(1566) Icarus and found da/dt = (—7. 3 ± 3.9) x 10~ 4 AU/Myr. Our best estimate is 



da/dt 



-3.2 ± 2.0) x 10~ 4 AU/Myr. 



Sitarskil (119981 ) found it necessary to incorporate 



a non-gravitational term da/dt = —58 x 10 4 AU/Myr in his orbit determination of 
(4179) Toutatis, however the availability of radar ranges in 1992, 1996, 2004, and 2008 



strong ly suggest a drift magnitude that does not exceed —5 x 10 4 AU/Myr. 



Ziolkowski 



(119831 ) examined the orbits of 10 asteroids and found drifts in four as t eroids , including 



a (-295.7 ± 14.6) x 10~ 4 AU/Myr drift for (1862) Apollo. 



Yeomand (119911 ) used a 



cometary model to search for perturbations and also detected a drift associated with (1862) 
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Apollo, though a value was not reported. Our best estimate is (—2.38 ± 0.25) x 10~ 4 
AU/Myr (Section [3]). It appears that these early estimates are not aligned with modern 
determi nations, and may hay e been caused by erroneous or insufficient astrometry. More 



recently, 



Chesley et al. 



(120081 ) searched for Yarkovsky signatures and reported rate estimates 



for 12 candidates. 



Here we use new developments in star catalog debiasing (I Chesley et al.ll2010l ) as well as 
the most recent astrometric data to compute semi-major drift rates for select NEAs, which 
multiplies the number of existing measurements by a factor of ~4. 

Observations of Yarkovsky rates can be used to place constraints on composition (i.e. 
metal vs. rock), physical properties (i.e. bulk density), and spin properties (i.e. prograde 
vs. retrograde). The magnitude of the force is dependent on the object's mass, size, 
obliquity, spin rate, and surface thermal properties. Separating how each of these quantities 
uniquely contributes to a measured da/dt is often not possible, but past Yarkovsky 
detections have allowed for insight into the associated objects. With certain assumptions 
on surface t hermal properties, b ulk densities were deter mined from the measured drifts 



of Golevka ( jChesley et al. 



20031 ) and (152563) 1992 BF flVokrouhlickv et al. 



2008). For 



the latte r, the magnitude and dir ection of the drift point to an obliquity in excess of 120 



degrees ( jVokrouhlicky et al. 



2008|). 



2. Methods 

2.1. Yarkovsky sensitivity 

The Yarkovsky drift manifests itself primarily as a change in mean anomaly (or 
along-track position), and some observational circumstances are poorly suited to detect 
such changes. Examples include optical astrometry secured when the line-of-sight is 
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roughly parallel to the asteroid velocity vector or when the object is at large distances from 
Earth. In both instances the differences in astrometric positions can be much smaller than 
observational uncertainties, resulting in low sensitivity to the Yarkovsky effect. The overall 
Yarkovsky sensitivity depends on the orbital geometry of the NEA and on the entire set 
of available observations. This can be quantified rigorously. For each epoch U at which 
optical observations were obtained (1 < % < N), we predict the position Pf for the best-fit 
orbit [dajdt = 0) and the position P* for the same orbit modified by a nominal non-zero 
da/dt. The value of the nominal rate is not important as long as it results in detectable 
(~arcsecond) changes in coordinates and as long as it is applied consistently to all objects; 
we used da/dt=0.1 AU/Myr. 

We then define the Yarkovsky sensitivity sy as 



where <7j is the positional uncertainty associated with observation i. This root mean square 
quantity provides an excellent metric to assess the relative sensitivity of any given data set 
to a drift in semi-major axis, including drifts caused by Yarkovsky influences. The metric 
can be applied to the entire set of available observations, or to the subset of observations 
that survive the outlier rejection steps described below. We computed both quantities and 
used the latter for our analysis. We found that data sets with scores sy below unity yield 
unreliable results, including artificially large rates and large error bars. Out of ~ 1,250 
numbered NEAs, only ~300 have sy > 1 and ~150 have sy > 2. In this paper we focus on 
a subset of these NEAs. 




N 



(3) 
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2.2. Orbital fits 



For this work we employed orbital fits to optical astrometry to determine semi-major 
axis drift rates for NEAs. We used the OrbFit software package , which is developed and 



maintained by the OrbFit Consortium (jMilani fc Gronchi 



20091 ). OrbFit can fit NEA 



trajectories to astrometric data by minimizing the root mean square of the weighted 
residuals to the data, optionally taking into account a given non-zero rate of change in 
semi-major a xis da/dt. We include d perturbations from 21 asteroids whose masses were 



estimated by 



Konopliv et al. 



(1201 lh . 



We downloaded optical astrometry for all numbered minor planets (NumObs.txt. gz) 
from the Minor Planet Center (MPC) on January 31st, 2012. We have assumed that 
all the astrometry has been properly converted to the J2000 system. The quality of the 
astrometry varies greatly, and we applied the data weighting and debias ing techniques 



Chesley et al. 



i mplem ented in OrbFit, which appear to follow the recommendations of 
( 120101 ) . Data weights are based on the time the observation was performed, the method 
of the observation (CCD or plate), the accuracy of the star catalog, and in some cases the 
accuracy of the observatory. Correction for known star catalog biases was applied when 
possible. Biases vary depending on the specific star catalog and region of the sky, and can 
reach 1.5 arcseconds in both right ascension (R.A.) and declination (Dec). Correction for 
these biases can substan tially improve t he rec overy of orbital parameters from observations. 



However, as discussed in 



Chesley et al. 



( 1201 Oh . not every observation can be debiased. Some 



observations were re ported to the MPC w ithout noting the star catalog used in the data 
reduction. Although 



Chesley et al. 



(120101 ) deduced the star catalogs used by several major 



surveys, there remain observations from smaller observatories that do not have associated 



star catalogs. Accordingly, a fraction of t 



Based on counts published 



l e astr ometry used in this paper was not debiased. 



Chesley et al. 



(120101 ). we estimate this fraction to be less than 
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7.2% of all the observations. 

Our procedure for determining the semi-major axis drift rate included three steps: an 
initial fit to the debiased data, an outlier rejection step, and a search for the best-fit da/dt, 
with iteration of the last two steps when necessary. 

We used the orbital elements from the Minor Planet Center's MPCORB database as 
initial conditions for the first fit for each object (step 1). This first fit, performed with 
da/dt = and outlier rejection turned off, slightly corrected the orbital elements for our 
weighted, debiased observations. The orbital elements from each object's first fit became 
the starting orbital elements for all later fits of that object. 

The second fit of each object served to reject outliers and was initially performed 
with da/dt = (step 2). The residual for each observation was calculated using the usual 
observed (O) minus computed (C) quantities: 

(R.A-o - R.A.c) x cos(Dec. ) \ 2 + / Deep - Dec. c \ 2 ^ 

where ctr.a. and o"Dec. are the uncertainties for that observation in R.A. and Dec, 
respectively. We rejected observations when their x res > a/8, and recovered previously 
rejected observations at Xrcs = \m , with the rejection step iterated to convergence. 
Results are fairly robust over a large range of thresholds for rejection (Section |3]). If the 
post-fit residuals were normally distributed, the chosen thresholds would result in < 1% of 
observations being rejected as outliers. Because errors are not normally distributed, our 
typical rejection rates are 2-5% of all available astrometry. This second fit produced the set 
of observations which were used in the third step. 

The third step was a series of orbital element fits to the observations over a set of fixed 
da/dt values. During these fits, we used the set of observations defined by the second fit 
and did not allow further outlier rejection. The quality of a fit was determined by summing 
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owest y 2 , we use d a 



19921 ) . A 



the squares of residuals x 2 — Yl Xrcs- To locate the region with the 
three-point parabolic fit or the golden-section minimization routine ( jPress et al.l 
parabola was then fit to the \ 2 curve in the vicinity of the minimum, and we used the 
minimum of the parabola to identify the best-fit da/dt value. 

Confidence limits were estimated using x 2 statistics. Confidence regions of 68.3% and 
95.4% (la and 2a, respectively) were established by the range of da/dt values that yielded 
X 2 values within 1.0 and 4.0 of the best-fit x 2 value, respectively (Fig. [Q. 

The initial outlier rejection step can in some cases eliminate valid observations simply 
because the Yarkovsky influences are not captured in a dynamical model with da/dt = 0. 
To circumvent this difficulty, we iterated the outlier rejection step with the best-fit da/dt 
value and we repeated the fitting process. In 52 out of 54 cases, the new best-fit value 
matched the previous best-fit value to within la, and we accepted the new best-fit values as 
final. For the other objects we repeated the reject and fit processes until successive best-fit 
values converged within la (which never required more than one additional iteration). Our 
results report the da/dt values obtained at the end of this iterative process. 



2.3. Sample selection 

We restricted our study to numbered NEAs with the best Yarkovsky sensitivity 
(Equation ED, specifically s Y > 2 (Fig. [2D- 

We also chose to focus on objects with non-zero da/dt values by using a signal-to-noise 
ratio (SNR) metric, defined as the ratio of the best-fit da/dt to its la uncertainty. We 
accepted all objects with SNR > 1 (Fig. [2}. 

Some asteroids have observations that precede the majority of the object's astrometry 
by several decades and have relatively high uncertainties. In order to test the robustness 
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of our results, we removed these sparse observations, which were defined as ten or fewer 
observations over a 10-year period. Fits were then repeated for these objects without the 
early observations. If the initial best-fit value fell within the la error bars of the new 
best-fit value, the initial result was accepted, otherwise, the object was rejected. 

Superior detections of the Yarkovsky effect are likely favored with longer observational 
arcs, larger number of observations, and good orbital coverage. For this reason we limited 
the sample to those NEAs with an observational arc at least 15 years long, with a number 
of reported observations exceeding 100, and with at least 8 observations per orbit on at 
least 5 separate orbits. 

We report on the 54 objects that met all of these criteria: sensitivity, SNR, sparse test, 
and orbital coverage. 

2.4. Validation 

We validated our optical-only technique whenever radar ranging observations were 
available on at least two apparitions. This could only be done for a fraction of the objects 
in our sample. In the remainder of this paper, optical-only results are clearly distinguished 
from radar + optical results. For the radar +optical fits, we included all available radar 
astrometry and disallowed rejection of potential radar outliers. The internal consistency of 
radar astrometry is so high that outliers are normally detected before measurements are 
reported. 

We also verified that a fitting procedure that holds successive da/dt values constant is 
equivalent to performing 7-parameter fits (6 orbital parameters and da/dt simultaneously). 
The da/dt values obtained with both procedures are consistent with one another. 
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2.5. Yarkovsky modeling 



In addition to the measurements described above, we produced numerical estimates of 
the diurnal Yarkovsky drift for each of the objects in our sample. Comparing the measured 
and estimated rates provides a way to test Yarkovsky models. In some instances, e.g., 
robust observations irreconcilable with accurate Yarkovsky modeling, it could also lead to 
the detection of other non-gravitational forces, such as cometary activity. Our numerical 
estimates were generated as follows. At e ach timestep, we computed the diurnal Yarkovsky 



acceleration according to equation (1) of 



Vokrouhlickv et al 



spherical body, with the physical parameters (lOpeil et al. 



fl2000f ). which assumes a 



20101 ) listed in Table [T] and an 



assumption of 0° or 180° obliquity. We assumed that the thermal conductivity did not have 
a temperature depend ence, but found that ad ding a temperature-dependent term according 



to the prescription of 



Hiitter fc Komlel (120081 ) (K = K + i^T 3 , with K x = 0.0076) did 



not change our predictions by more than 1%. We then resolved the acceleratio n along 



ortho gonal directions, and used Gauss' form of Lagrange's planetary equations ( iDanby 
19921 ) to evaluate an orbit-averaged da/dt. 



The physical parameters chosen for these predictions mimic two extremes of rocky 
asteroids; one is intended to simulate a rubble pile with low bulk density, the other a 
regolith-free chunk of rock (Tabled]). These parameters cor respond to a therm al inertia 



Delbo et al. 



(120071 ). who found 



range of 77 — 707 J m -2 s -a5 K _1 , enveloping the results of 
an average NEA thermal inertia to be 200 J m -2 s -0 ' 5 K _1 . In most cases, the drift rates 
produced by these two extreme cases encompass the drift produced by a rubble-pile object 
that has a regolith-free surface, or the drift produced by a solid object with regolith. 

There is no simple relationship between these physical parameters and predicted drift 
rates, but for most cases the rubble pile exhibits the larger da/dt values due to its low bulk 
density (Equation [2]). The smaller values of density of the surface and thermal conductivity 
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Table 1: Physical and thermal properties used for numerical estimates of the semi-major 
axis drift of asteroids. Thermal propertie s are based on measurements of three meteorites at 



200 K, as measured by 



Opeil et al. 



( 120101 ). Listed are heat capacity C, thermal conductivity 



K, bulk density of the surface p s , and mean bulk density p. 



Composition 


C (J kg^K- 1 ) 


K (W rn^K- 1 ) 


Ps (kg m 3 ) 


P (kg m 3 ) 


Rubble Pile 


500 


0.01 


1200 


1200 


Rock Chunk 


500 


0.50 


2000 


2000 



for rubble piles produce a smaller thermal inertia, and therefore a longer thermal lag. 
Generally, but not always, this longer thermal lag, combined with the rotation of the 
asteroid, allows for a larger fraction of departing thermal emission to be aligned with the 
asteroid's velocity, resulting in a larger drift. 

When available, measured values of the geom etric albedo, diameter, and spin rate from 



the JPL Small-Body Database (jChamberlin 



20081 ) were incorporated into our predictions 



for Yarkovsky drifts. When not available, the diam eter D in km was estimated from the 



absolute magnitude H using ( iPravec h Harris 



20071), 



D 



1329 



* 10 



-0.2H 



(5) 



where we used two values of the V-band geometric albedo pv (0.05 and 0.45), a range that 
captures observed albedos for the majority of NEAs. When spin rate was unknown, we 
assumed a value of 5 revolutions/day, ba sed on the average spin rate values for asteroids 1 
to 10 km in diameter shown in Fig. 1 of IPravec fc Harris! ( 120001 ) . Emissivity was assumed 
to be 0.9. Bond albedo was estimated with a uniform value of the pha se integral (q=0.39 



on the basis of the IAU two-parameter magnitude system for asteroids 
and an assumed slope parameter G=0.15. 



Bowell et al. 



(1989) 
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We have assumed pv = 0.14 for the purpose of quantifying the Yarkovsky efficiency 
when the asteroid size was unknown. 

3. Results 

We measured the semi-major axis drift rate of all 1,252 numbered NEAs known as of 
March 2012. Some of the drift rates are not reliable because of poor sensitivity to Yarkovsky 
influences (Fig. [2]). 

After our process of selection and elimination (Section I2.3[) . we were left with 54 NEAs 
that exhibit some of the most reliable and strongest drift rates. Although we report objects 
with sy > 2, we have the most confidence in objects with highest Yarkovsky sensitivity, 
and we show objects in order of decreasing sy value in our figures. 

We examined the impact of various choices of reject /recover thresholds when rejecting 
outlier observations (Fig. [3]). At moderate values of the rejection threshold (i.e. eliminating 
less than ~5% of observations), best-fit values are consistent with one another. In this 
regime, results are fairly robust against the choice of rejection thresholds. However results 
do become sensitive to rejection thresholds when a larger fraction of observations is rejected. 
As the reject /recover thresholds become more stringent, astrometry with evidence of 
semi-major axis drift is preferentially rejected, and the best-fit da/dt values approach zero. 
Our adopted reject /recover thresholds (y/8/y/7) are stringent enough that they eliminate 
obvious outliers, but not so stringent as to suppress the Yarkovsky signal. In 52 out of 
54 cases, repeating the outlier rejection step with the best-fit da/dt value resulted in no 
appreciable change to the result. 

As a validation step, we compared the semi-major axis drift rates obtained with our 
procedure (both optical-only and radar +optical) to previously published values (Table EJ). 
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We found good agreement for Go 



1992 BF flVokrouhlickv et al. 



study done by 



Chesley et al. 



Chcs ley et al 



(12008 



evka flCheslev et all 120031 : ICheslev et al.l 12008( 1 and 
2008), and for most, but not all, NEAs included in a similar 
(120081 ). The differences between our results and those of 
) can probably be attributed to our use of debiased data, of improved 
data weights, and of longer observational arcs extending to 2012. Eight objects included 
in Table [2] meet our selection criteria for detailed analysis in the rest of this paper: 
(1620) Geographos, (1685) Toro, (1862) Apollo, (1865) Cerberus, (2063) Bacchus, (2100) 
Ra-Shalom, (2340) Hathor, and (152563) 1992 BF. 

Several conclusions can be drawn from the data presented in Table [2J First, the RMS 
values indicate excellent fits to the astrometry. Second, the solutions with non-zero da/dt 
values provide a much better match to the radar data than the gravity-only solutions, with 
typical RMS values decreasing by a factor of 2 or more. Third, radar -l-optical estimates 
have consistently lower error bars than optical-only estimates, sometimes dramatically so, 
which is typical in NEA studies. Finally, there is a generally good agreement between 
the optical-only da/dt values and the radar -l-optical da/dt values, indicating that the 
optical-only technique is a useful tool that can be used even in the absence of radar data. 

Drift rates for the 54 NEAs that pass our selection criteria are presented in Table |3] 
along with orbital elements and physical properties. If an object has both a optical-only 
and a radar -l-optical value, we used the more accurate radar -l-optical value in the following 
figures and calculations (unless specified otherwise). We used Equation (T5]) with a density 
of 1,200 kg m -3 to compute efficiency factors fy and found that objects divided roughly 
into two groups. 

In the first group of 42 objects with fy < 2 x 10 -5 , most observed da/dt values are 
consistent (within la) with Yarkovsky predictions. We refer to these objects as Yarkovsky- 
dominated (Figs. @]and|5]). In the second group of 12 objects with fy > 2 x 10 -5 the 
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observed da/dt values are somewhat larger than Yarkovsky predictions, but improvements 
in the knowledge of physical properties or in Yarkovsky modeling could plausibly bring 
some of the observed rates in agreement with predictions. We refer to these objects as 
possibly Yarkovsky- dominated (Fig. [6]). 

Figures H] and indicate that there is generally agreement between observations and 
numerical estimates of Yarkovsky drift rates for NEAs with /y<2x 10~ 5 . These data 
suggest that fy ~ 10~ 5 represents a typical efficiency for the Yarkovsky process. Predicted 
values are based on calculations with obliquities of 0° and 180°, therefore, observed rates 
that are lower than predictions could still be due to the Yarkovsky effect. 

The majority of objects in Fig. [5] appear to exceed predictions. This is a consequence 
of the SNR > 1 selection criterion, as it eliminates objects with lower da/dt values. 

On the basis of the entire sample of measured drifts for objects with sy > 2, we can 
compute average properties for observed Yarkovsky rates and efficiencies. The mean, mean 
weighted by uncertainties, median, and dispersion are shown in Table HI The aggregate 
properties are comparable if we restrict objects to the subset with SNR > 1, except for 
slightly increased da/dt rates (median rate of ~ 12 x 10~ 4 AU/Myr instead of ~ 7 x 10~ 4 
AU/Myr), as expected. The Yarkovsky process appears to have an efficiency fy of order 
10~ 5 , with a fairly small dispersion. Because the Yarkovsky efficiency scales with density 
(fy\ P = /y|i,2oo x 200 kg m~ 3 ) some of the observed scatter is due to density variations. 

4. Discussion 

In this section we examine several consequences of our results. First we discuss how 
the Yarkovsky drifts can inform us about asteroid physical properties, spin states, and 
trajectories. Then we discuss binary asteroid (1862) Apollo and the curious case of asteroid 
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(1036) Ganymed. Finally we discuss the possible mechanisms for non-Yarkovsky driven 
rates, including association with meteoroid streams and rock comet phenomenon. 



4.1. Yarkovsky-derived constraints on asteroid physical properties 



Because a clear connection exists between asteroid physical properties and Yarkovsky 
drifts, we explored the constraints that can be placed on bulk density and surface thermal 
conductivity for seven objects with well-known diameters and (excepting one case) spin 
periods: (1620) Geographos, (1862) Apollo, (2100) Ra-Shalom, (2062) Aten, (2340) Hathor, 
(1566) Icarus, and (3361) Orpheus. We compared the measured Yarkovsky rates to 
numerical estimates obtained with a range of physical parameters. For these estimates, 
we assumed a constant heat capacity C = 500 J kg~ x K _1 (Tabled]) and a single value of 
the bulk density of the surface p s — 1, 700 kg m -3 , but we explore a wide range of bulk 
density and surface thermal conductivity values. Because the obliquities are uncertain or 
ambiguous in many cases, we chose to illustrate outcomes for two obliquity values, typically 
180° and 135°. 



Our results are shown in Figures [7] and El which are similar to Fig. 4 in IChesley et al 



(120031 ). The shaded range consistent with the la confidence limits on da/dt delineates the 
space of acceptable bulk densities and thermal conductivities, assuming that the Yarkovsky 
effect is being modeled correctly. By acceptable, we mean consistent with observed da/dt 
values, even though some of the K — p values may not be appropriate for asteroids. 



Infrared observations indica te that (2100) Ra - Shalom has a 



between 0.1 and 1 W m K (iDelbo et al 



2003 



Shepard et al 



therm al conductivity 



20081 ). If we assume a 



minimum bulk density of 1,500 kg m 3 , this conductivity value is consistent with the range 
suggested by our Yarkovsky rate determination. 
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If we make the same minimum density assumption for (1620) Geographos, our 
measurements suggest that its surface thermal conductivity is greater than 0.002 
W m- 1 K" 1 . 

For (1862) Apollo, we show the range of physical properties that are consistent 
with both the optical-only fits and the radar+optical fits. The precision of the radar 
measurements dramatically shrinks the size of the measured error bars, with correspondingly 
tighter constraints on density and surface thermal conductivity. This example illustrates 
that reliable obliquity determinations will be important to extract physical properties from 
Yarkovsky rate determinations. 

Our measurement of (2062) Aten's drift provides some useful insights. If we assume 
that its bulk density exceeds 1,500 kg m~ 3 , then its surface thermal conductivity K must 
exceed 0.3 W var 1 K _1 . Furthermore, if we assume that its bulk density exceeds 1,600 
kg m~ 3 , the la confidence region on the measured Yarkovsky drift suggests that its obliquity 
is between 180° - 135°. 

The Yarkovsky simulations for (2340) Hathor were computed with an assumed spin 
period of 4.5 hours. If the actual period is longer, the curves shown would shift to the 
left, and if the period is shorter, the curves would shift to the right. Consequently, we 
cannot make inferences about the K value for this object until its spin period is measured. 
However, looking at the height of the curve, and with an assumption that the object's bulk 
density is greater than 1,500 kg m~ 3 , we can conclude that (2340) Hathor likely has an 
obliquity lower than 180°. 

The assumption of 135° or 180° obliquity for (1566) Icarus restricts this object to low 
surface conductivity values and low bulk density values, or high surface conductivity values 
and high bulk density values. Although these obliquities do produce physically plausible 
parameter combinations, it seems likely that the obliquity for this object is < 135°. 
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The curves for (3361) Orpheus were calculated with an assumed geometric albedo 
of 0.15. As (3361) Orpheus has a positive da/dt value, obliquities were assumed to be 
0° and 45°. The curve representing an obliquity equal to 0° for this object requires very 
low (< 0.002 W m _1 K _1 ) or very high (> 0.7 W m _1 K _1 ) surface thermal conductivity 
values for most densities. A more likely scenario is that this object has an obliquity > 0°, 
or perhaps even > 45°. An independent measurement of the obliquity could be used to 
validate obliquity constraints derived from Yarkovsky measurements. 



4.2. Yarkovsky rates and distribution of spin states 



La Spina et al.1 (120041 ) and IChesley et al.l (120081 ) examined the predominance of 
retrograde spins and negative Yarkovsky drift rates and concluded that they were consistent 
with the presumed delivery method of NEAs from the main belt of asteroids . The z/ 6 and 



3:1 resonance regions deliver NEAs to near-Earth space (iBottke et al.l l2002a[ ). A main belt 
asteroid can arrive at the 3:1 resonance at 2.5 AU via a positive (if it originates in the inner 
main belt) or negative (if it originates in the outer main belt) Yarkovsky drift. However, a 



main belt asteroid can only arrive at t he ve, resonan c e (at t 



way o f a negative drift. According to 



Bottke et al. 



re in ner edge of the main belt) by 



(I2002af ) and 



Morbidelli fc Vokrouhfickv 



( 120031 ). 30% — 37% of NEAs are transported via the uq resonance, with the rest from other 
resonances. The net result is a preference for retrograde spins. 

An observational consequence of th is process would be an excess of retrograde rotators 
in the near-Earth asteroid population. lLa Spina et al.l (120041 ) conducted a survey of 21 
NEAs and found the ratio of retrograde/prograde rotators to be 2.0lg 7 . 

We note that out of the 42 Yarkovsky-dominated NEAs, 12 have a positive da/dt value. 
For this sample, our ratio of retrograde/prograde rotators is 2.5 ± 0.1, similar to the value 
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found by lLa Spina et al.l ( 120041 ) 



4.3. Impact of drift rates on asteroid trajectory predictions 

The semi-major axis drifts described in this paper affect NEA trajectory predictions. 
An orde r of magnitude estimate fo r the along track displacement due to a non-zero da / dt is 



given m 



Vokrouhlicky et al. 



feOQOh : 



Ap ~ 7a 4 (A 10 t) 2 aJJ 2 



(6) 



where Ap is in units of km, d 4 is da/dt in 10 -4 AU/Myr, A 10 t is the time difference 
between observations in tens of years, and ciau is the semimajor axis of the object in 
AU. For instance, the estimated along-track displacement due to the observed da/dt for 
(1862) Apollo is 9 km after 10 years. Similarly, the estimated along-track displacement for 
faster-moving (1864) Daedalus is 67 km after 10 years. 

Our data indicate that (101955) 1999 RQ36, the target of the OSIRIS-REx mission, has 
a measurable Yarkovsky drift of (—18.9 ± 0.2) x 10~ 4 AU/Myr. Although it has a relatively 
short arc (12 years) it has been observed three times by radar, allowing for an accurate 
da/dt measurement. We estimated the along-track displacement of (101955) 1999 RQ36 
over the 6-month duration of the OSIRIS-REx mission to be 0.3 km, which will be easily 
detectable by a radio science instrument. 



4.4. Binary asteroid (1862) Apollo 



(1862) Apollo is a binary asteroid (jOstro et al.ll2005l ). Binary asteroids present a 
unique opportunity for the determination of physical parameters. If mass and density 
can be measured from the binary orbit and component sizes, the Yarkovsky constraint on 
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thermal conductivity can become much more meaningful. If the orientation of the plane of 
the mutual orbital can be measured, a plausible obliquity can be assumed, which makes the 
constraints on thermal properties tighter still. In some cases, actual obliquity measurements 
can be obtained from shape modeling efforts. 



Yeomansl (11991 



19921 ) identified a non-gravitational perturbation acting on the orbit of 
(1862) Apollo, but was not able to determine a drift magnitude. To la our observed da/dt 
value for (1862) Apollo agrees with our Yarkovsky predictions. 



4.5. The curious case of (1036) Ganymed 

(1036) Ganymed has by far the largest Yarkovsky efficiency value (/y ~ 15 x 10 -5 ) 
among the objects presented in Table [31 With a nominal value of ~ —7 x 10 -4 AU/Myr, 
the measured da/dt value is comparable to that of other NEAs. Combined with Ganymed's 
large diameter estimate (~ 32 km based on IRAF measurements), this Yarkovsky rate 
results in an unusually high f Y value. 

How can this anomaly be explained? One possibility is that some of the early 
astrometry, dating back to 1924, is erroneous. This could be due to measurement errors, 
timing errors, bias errors, or reference frame conversion errors. We evaluated the semi-major 
axis drift with various subsets of the available astrometry and found values ranging between 
—3 x 10~ 4 and —8 x 10 -4 AU/Myr. On that basis we modified the adopted uncertainties 
for this object, and our preferred value is (—6. 62^^) x 10 -4 AU/Myr. Doing so does not 
eliminate the possibility of systematic bias in the astrometry, and we are still left with 
anomalously high /y values. 

Another possibility is that the diameter of Ganymed, an S-type asteroid, is much 
smaller than reported. This seems unlikely considering the more recent WISE albedo 
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measurement of pv = 0.212 (IMasiero et al.ll2012l ) which suggests a diameter of ~ 36 km. 



If Ganymed's bulk density was especially low, a higher than usual fy value would be 
expected, but this would likely explain a factor of 2 or 3 at most, and would not explain the 
anomalous value. 

Perhaps Ganymed departs significantly from a spherical shape, with an effective 
diameter and mass that are much smaller than those implied by the diameter values 
reported in the literature. The relatively low lightcurve amplitudes do not seem to support 
such an argument, unless the asteroid is particularly oblate. In that case one could plausibly 
arrive at volume and mass estimates that are off by a factor of 5-10. 

If we can rule out these possibilities (i.e. Ganymed is roughly spherical with no 
substantial concavities, its diameter estimate is reasonably accurate, and the early 
astrometry can be trusted), and if no other modeling error can be identified, then we would 
be compelled to accept an anomalously high Yarkovsky efficiency for this object. 



4.6. Non- Yarkovsky processes 

In the course of our study we observed drift values that cannot be accounted for 
easily by Yarkovsky drift, because they considerably exceed the predicted Yarkovsky 
rates. In most cases, these can be attributed to poor sensitivity to Yarkovsky influences 
(Fig. [2]). Therefore, the high rates can generally be safely discarded. In other cases, the 
high rates may be due to erroneous optical astrometry or mismodeling of asteroid-asteroid 
perturbations. However we cannot entirely rule out the possibility that some of the high 
drift rates are secure and will be confirmed by further observation and analysis. If the high 
rates cannot be ascribed to poor Yarkovsky sensitivity or faulty astrometry, one would need 
to invoke other non-gravitational forces. 
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One possibility is that orbits are perturbed when NEAs are losing gas or dust in an 
anisotropic manner. To estimate a rough rate of mass loss that would be needed to account 
for the drifts measured, we used the basic thrust equation 



F = qV e 



(7) 



where F is the force, q is the rate at which the mass departs the asteroid, and V e is the 
ejection speed. For an asteroid of mass m this yields 



^mass loss 



qVe 
m 



which can be incorporated into Gauss' form of Lagrange's planetary equations (IDanby 
19921 ) as an acceleration aligned with the velocity of the object. The dependence of the 
force on heliocentric distance r is not known precisely; we assumed F oc r -2 , similar to the 



Yarkovsky dependence, for simplicity, and because the amount o 



with the amount of incident radiatio n (as in Fig. 4 of 



V e — 1.5 m s 1 , the value derived by 



HsiehetaL 



outga ssing likely scales 



Delsemmd (119821 ) ) . We assumed 



(120041 ) for 133P/Elst-Pizarro, and we 



assumed that the mass is departing in the optimal thrust direction. 

We quantified the mass loss rates needed to produce the observed drifts of NEAs with 
the highest Yarkovsky efficiencies. We estimated a rate of 0.16 kg s _1 for (154330) 2002 
VX94 and 2.3 kg s^ 1 for (7889) 1994 LX. Although these estimates represent the minimum 
amount of mass loss necessary to account for the observed drifts (if due to mass loss), they 
are smaller than typical levels from comets. Comets have mass loss rates that span a wide 



range of values. On 



( jJewitt fc Matthews 



;he h igh side a rate of 2 x 10 6 kg s 1 was est imated for Hale-Bopp 



19991). On the low side 



Ishiguro et al. 



(120071 ) measured mass loss 



rates for three comets, averaged over their orbits: 2P/Encke (48 ± 20 kg s 1 ), 22P/Kopff 
(17 ± 3 kg s" 1 ), and 65P/Gunn (27 ± 9 kg s^ 1 ). Mass loss rates of active asteroids have 



been estimated to be in the range 



(107P/Wilson-Harrington) (jJewitt 



rom < 0.04 kg s^ 1 (113P/Elst-Pizarro) to < 150 kg s" 1 



20121 ) . 
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Mass loss does not seem to be a viable mechanism to explain the semi-major axis drift 
rate of (1036) Ganymed, as it would require a minimum mass loss rate of ~2,500 kg s _1 . 
This would presumably have left detectable observational signatures, which have not been 
reported to date. 

We explore a couple of possibilities for mass loss mechanisms that could cause 
semi- major axis drifts. 



4-6.1. Associations with meteoroid streams 



To our knowledge, (433) Eros, (1566) Icarus, (1620) Geographos, (1685) Toro, (1862) 



Apollo, and 1982 ' 
meteoroid stream. 



'A are th e only 



objects in our sample to have been associated with a 



Sekaninal (119761 ) found a weak correlation between the first five objects 



and various streams using the 
described as not convincing by 



dissimilari t y crit erion" . However, this metric was later 



Jenniskend (120081 ). and current literature does not support 



such associations. In our results, Apollo shows good agreement with Yarkovsky predictions, 
with /y = 0.25 x 10 -5 . The Yarkovsky force is therefore a plausible cause of Apollo's 
observed semi- major axis drift. 



4-6.2. Rock comet phenomenon 



The brightening of (3200) Phaethon, the parent body of the Gem inid meteor shower 



has been attributed to a "rock comet" phenomenon (IJewitt fc Li 



20101 ) . Wit h a perihelion 



Jewitt fc Li 



at 0.1 4 AU, (3200) Phaethon's surface temperatures have been estimated by 
(120101 ) to be in the range 746 < T < 1050 K. The authors propose that these high surface 
temperatures could create thermal gradients in the body, resulting in thermal fracturing 



that would release dust. The resulting mass loss would affect the orbit. The combination 
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of mass loss due to decomposing hydrated minerals and thermal fracturing led the authors 
to term (3200) Phaethon a "rock comet". A moderate amount (~ 1 kg s -1 ) of mass lost in 
an anisotropic manner by "rock comets" could explain the observed semi-major axis drift 
rates. 



5. Conclusions 

Modeling of the Yarkovsky effect is needed to improve trajectory predictions of 
near-Earth asteroids and to refine our understanding of the dynamics of small bodies. 
Using fits to astrometric data, we identified semi-major axis drifts in 54 NEAs, 42 of which 
show good agreement with numerical estimates of Yarkovsky drifts, indicating that they 
are likely Yarkovsky-dominated. These objects exhibit Yarkovsky efficiencies of ~10 -5 , 
where the efficiency describes the ratio of the change in orbital energy to incident solar 
radiation energy. 12 objects in our sample have drifts that exceed nominal Yarkovsky 
predictions and are labeled possibly Yarkovsky-dominated. Improvements in the knowledge 
of physical properties or in thermal modeling could bring these drift rates in better 
agreement with results from numerical models. However, if the high rates are confirmed 
by additional observations and analysis, they would be indicative of the presence of other 
non-gravitational forces, such as that resulting from a loss of mass. 

None of this work would have been possible without the availability of the OrbFit 
software package (available at http://adams.dm.unipi.it/orbfit/). 

We are grateful to The Minor Planet Center and all astronomers who submitted data 
to the Minor Planet Center. 

We thank Ned Wright (UCLA) for insightful comments. 
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-0.0005 
da/dt (AU/myr) 



0.0005 



Fig. 1. — Search for best-fit da/dt value to optical astrometry of (2100) Ra-Shalom (1,281 
observation epochs, 2,562 observations, 7 adjustable parameters, 2,555 degrees of freedom). 
The sums of squares of residuals corresponding to a range of da/dt values are shown as circles, 
with a parabolic fit shown as a dotted li ne. The da/dt valu es plotted here were determined 



by the golden section search algorithm (jPress et al. 



1992T ) as it searched for and found a 



minimum at da/dt = —5.20 x 10 4 AU/Myr with a reduced \ 2 value of 0.30. Confidence 
limits of 68.3% (la) are indicated by the thick dashed line, and correspond to the range 
da/dt = [-7.4, -2.9] x 10~ 4 AU/Myr. The thin dashed line shows the 95.4% (2a) confidence 



region. 
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-0.02 0.02 

Semi — major axis drift (AU/Myr) 



Fig. 2. — Yarkovsky sensitivity metric sy plotted as a function of semi-major axis drift rate 
da/dt for 1,252 numbered NEAs. Data sets with Yarkovsky sensitivity below unity (dashed 
line) yield unreliable results, including large rates and large error bars. Our selection criteria 
require sy > 2 (dotted line) and SNR > 1. The 80 objects that meet both selection criteria 
are shown in green. About 26 of these 80 NEAs are eliminated by the sparse test and orbital 
coverage requirements (see Section l2~3l) . 
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Fig. 3. — Impact of different choices of reject /recover thresholds for the initial rejection step 
(da/dt = 0) on the best-fit da/dt values. Results from optical-only fits are shown with their 
la error bars for two representative cases, (2202) Pele and (2063) Bacchus. Best-fit da/dt 
values are consistent with one another in the left half of the diagram. Values to the right 
of 2.3/2.2 (Pele) and 2.5/2.4 (Bacchus) have SNR less than unity and would not meet our 
selection criteria. Our adopted reject /recover thresholds are \/8 = 2.828 and y/l = 2.646. 
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Table 2: Comparison of our o ptical-only resu 



results of a previous study by 



Chesley et al- 



ts to radar +optical (r+o) results and to the 
(j2008l ). Best-fit da/dt values in units of 10 -4 



AU/Myr and their one-sigma uncertainties are listed for optical-only and radar -l-optical ob- 
servations. Also shown is the root-mean-square (RMS) of weighted residuals for the gravity- 
only (da/dt = 0) solution and for the non-zero da/dt solution (RMS'). We restrict the radar 
analysis to those objects that have range measurements on at least two apparitions; this 
excludes (1685) Toro, (1865) Cerberus, (2063) Bacchus, (2340) Hathor, (85953) 1999 FK21, 
and (152563) 1992 BF. 



NEA 



radar radar r+o 
RMS RMS' RMS' 



optical-only 

da/dt 



r+o Chesley 08 

da/dt da/dt 



1.111 0.403 



(1620) Geographos 0.393 0.356 

(1685) Toro 

(1862) Apollo 
(1865) Cerberus 
(2063) Bacchus 
(2100) Ra-Shalom 
(2340) Hathor 
(6489) Golevka 
(54509) YORP a 

(85953) 1999 FK21 b 

(101955) 1999 RQ36 b 15.694 0.127 
(152563) 1992 BF C 



0.488 0.594 

0.879 0.387 
0.796 0.260 



0.55 
0.51 
0.61 
0.54 
0.59 
0.51 
0.67 
0.61 
0.55 
0.56 
0.39 
0.60 



-2.43 + 0.7 
-1.40 + 0.7 
-1.79 + 0.6 
-5.11 + 2.7 
-4.17 + 3.7 
-4.79 + 2.2 
-14.55 + 3.6 
-2.05 ± 12.6 
-25.98 + 37.4 
-10.44+1.5 
-12.90 ± 7.1 
-1 2.81 ± 1.0 



-2.52 + 0.6 



-2.38 + 0.3 



-5.45 + 1.5 

-5.74 + 0.7 
-35.63 + 10.5 

-18.9 + 0.2 



-1.18 + 0.39 
-0.52 + 0.27 
-2.44 + 0.26 
-7.80 + 2.28 
-10.59 + 2.21 
-7.09 + 0.88 
-13.94 + 3.84 
-6.39 + 0.44 
-25.12 + 6.18 
-14.13 + 2.35 
-15.69 + 4.99 
-10.78 + 0.73 



This object is in a Sun-Earth horseshoe orbit (ITaylor et al.ll2007). b This ob 



experiences perihelion precession of ~16 arcseconds /c entury (jMargot fc Giorgini 



This object is the target of the OSIRIS-REx mi ssion ( Chesle y et a 



object use the astrometry corrections given in 



Vokrouhlicky et al 



ect 



2010). 



2012 ). d Fits to this 



d2008|) for the 1953 



observations, which we did not subject to rejection. 
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Fig. 4. — Measured and predicted drift values for 20 asteroids with Yarkovsky-dominated 
drifts, ordered by decreasing value of Yarkovsky sensitivity. Best fits to optical-only data 
are shown as squares with dotted la error bars. Shaded boxes show a range of predicted 
Yarkovsky rates representing different compositions (Tabled]). As predicted Yarkovsky values 
were calculated assuming 0° or 180° obliquity, the shaded boxes represent maximum drifts 
for the object. Therefore, a fit that lies between a shaded box and da/dt = is considered 
to have good agreement. Objects with a single corresponding shaded box have a known 
diameter (Table [3]). Objects with two shaded boxes did not have known diameters, and were 
modeled using diameters derived from assumed albedos (45% in light blue, larger predicted 
drift magnitudes, and 5% in dark blue, smaller predicted drift magnitudes). The vertical 
extents of the shaded boxes represent the range of compositional types described in Table 
[TJ with the larger absolute values representing the "rubble pile" composition, and the lower 
absolute values representing the "rock chunk" composition. 
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Fig. 5. — Measured and predicted drift values for an additional 22 asteroids with Yarkovsky- 
dominated drifts, ordered by decreasing value of Yarkovsky sensitivity. Symbols are as in 
Fig. HJ The observed rates for the majority of objects shown in this figure appear to exceed 
predicted values. This is a consequence of the SNR > 1 selection criterion which eliminates 
objects with lower da/dt values. 



-32 - 



80 



— ] 1 r 



60 



■e- 40 



ZD 

i 20 



□ 



-20 




-40 



■o 

CD 

E 

CO 



CO 

m 



< 
I- 

oj 

00 



CO 



2- CO 



o> 
X 

> 

CM 
O 
O 
CNJ 



— ' 



CL 
CD 
N 



CO 
CD 

m 
o 



eg 

CO 
CO 



Fig. 6. — Measured and predicted drift values for 12 asteroids with possible Yarkovsky- 
dominated drifts, defined as objects with Yarkovsky efficiency /y exceeding 2 xl0~ 5 . Sym- 
bols are as in Fig. HI Most objects in this figure have measured drifts that lie outside of the 
range of values expected on the basis of Yarkovsky models. This could be due to inaccuracies 
in our knowledge of physical properties, faulty astrometry, or modeling errors. 
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Fig. 7. — Range of bulk densities and thermal conductivities of three Yarkovsky-dominated 
asteroids consistent with their observed da/dt values. Blue (top) solid line corresponds to 
values consistent with best-fit da/dt and 180° obliquity, pink (lower) solid line corresponds 
to values consistent with best-fit da/dt and 135° obliquity. Dashed regions surrounding each 
solid line encompass the la confidence limits on the corresponding da/dt determinations. 
Not all values displayed in this K-p space are necessarily appropriate for asteroids. Infrared 
observations suggest that (2100) Ra-Shalom has a thermal conductivity between 0.1 and 1 



W m- 1 K- 1 flDelbo et al. 



2003 



Shepard et al 



20081 ) , consistent with the range suggested by 



our Yarkovsky rate determination. For Apollo, we show results for both optical-only and 
radar+optical determinations. The inclusion of radar data greatly reduces the error bars on 
the measured drift, and therefore the area of the shaded curves. 
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(2062) Aten: optical-only (2340) Hathor: optical-only 




Surface thermal conductivity, K (W m" 1 K" 1 ) Surface thermal conductivity, K (W m" 1 K" 1 ) 



Fig. 8. — Companion to Fig. [7J Range of bulk densities and thermal conductivities of three 
Yarkovsky-dominated asteroids consistent with their observed da/dt values. For (2062) Aten 
and (2340) Hathor, blue (top) solid line corresponds to values consistent with best-fit da/dt 
and 180° obliquity, pink (lower) solid line corresponds to values consistent with best-fit 
da/dt and 135° obliquity. The constraints for (1566) Icarus suggest that it may have a lower 
obliquity than those assumed. (3361) Orpheus has a positive drift, so the blue (top) solid 
line corresponds to values consistent with best-fit da/dt and 0° obliquity, pink (lower) solid 
line corresponds to values consistent with best-fit da/dt and 45° obliquity. Dashed regions 
surrounding each solid line encompass the la confidence limits on the corresponding da/dt 
determinations. Not all values displayed in this K-p space are necessarily appropriate for 
asteroids. A period of 4.5 hours was assumed for (2340) Hathor, and a 0.15 geometric albedo 
was assumed for (3361) Orpheus. 



-35 - 
REFERENCES 

Bottke, W. F., Morbidelli, A., Jedicke, R., Petit, J., Levison, H. F., Michel, P., & Metcalfe, 
T. S. 2002a, Icarus, 156, 399 

Bottke, Jr., W. F., Vokrouhlicky, D., Rubincam, D. P., & Broz, M. 2002b, Asteroids III, 
395 

Bottke, Jr., W. F., Vokrouhlicky, D., Rubincam, D. P., & Nesvorny, D. 2006, Annual 
Review of Earth and Planetary Sciences, 34, 157 

Bowell, E., Hapke, B., Domingue, D., Lumme, K., Peltoniemi, J., & Harris, A. W. 1989, in 
Asteroids II, 524-556 

Chamberlin, A. B. 2008, Bulletin of the American Astronomical Society, 40 

Chesley, S., Ostro, S., Vokrouhlicky, D., Capek, D., Giorgini, J., Nolan, M., Margot, J. L., 
Hine, A., Benner, L., & Chamberlin, A. 2003, Science, 302, 1739 

Chesley, S. R. 2006, Asteroids, Comets, and Meteors conference abstracts 

Chesley, S. R., Baer, J., & Monet, D. G. 2010, Icarus, 210, 158 

Chesley, S. R., Vokrouhlicky, D., Ostro, S. J., Benner, L. A. M., Margot, J. L., Matson, 
R. L., Nolan, M. C, & Shepard, M. K. 2008, LPI Contributions, 1405, 8330 

Chesley, S. R. et al. 2012, Asteroids, Comets, Meteors conference abstracts 

Danby, J. M. A. 1992, Fundamentals of Celestial Mechanics (Willmann-Bell, Inc) 

Delbo, M., dell'Oro, A., Harris, A. W., Mottola, S., & Mueller, M. 2007, Icarus, 190, 236 

Delbo, M., Harris, A. W., Binzel, R. P., Pravec, P., & Davies, J. K. 2003, Icarus, 166, 116 

Delsemme, A. H. 1982, in Comets, ed. L. L. Wilkening, 85-130 



-36 - 

Giorgini, J. D., Benner, L. A. M., Ostro, S. J., Nolan, M. C, & Busch, M. W. 2008, Icarus, 
193, 1 

Giorgini, J. D., Ostro, S. J., Benner, L. A. M., Chodas, P. W., Chesley, S. R., Hudson, 
R. S., Nolan, M. C, Klemola, A. R., Standish, E. M., Jurgens, R. F., Rose, R., 
Chamberlin, A. B., Yeomans, D. K., & Margot, J. L. 2002, Science, 296, 132 

Goldreich, P. & Sari, R. 2009, The Astrophysical Journal, 691, 54 

Hsieh, H. H., Jewitt, D. C, & Fernandez, Y. R. 2004, The Astronomical Journal, 127, 2997 

Hiitter, E. S. & Komle, N. I. 2008, in 5th European Thermal-Sciences Conference, The 
Netherlands 

Ishiguro, M., Sarugaku, Y., Ueno, M., Miura, N., Usui, F., Chun, M.-Y., & Kwon, S. M. 
2007, Icarus, 189, 169 

Jenniskens, P. 2008, Earth Moon and Planets, 102, 505 

Jewitt, D. 2012, The Astronomical Journal, 143, 66 

Jewitt, D. & Li, J. 2010, The Astronomical Journal, 140, 1519 

Jewitt, D. & Matthews, H. 1999, The Astronomical Journal, 117, 1056 

Konopliv, A. S., Asmar, S. W., Folkner, W. M., Karatekin, O., Nunes, D. C, Smrekar, 
S. E., Yoder, C. F., & Zuber, M. T. 2011, Icarus, 211, 401 

La Spina, A., Paolicchi, P., Kryszczynska, A., & Pravec, P. 2004, Nature, 428, 400 

Margot, J. L. & Giorgini, J. D. 2010, in IAU Symposium, ed. S. A. Klioner, P. K. Seidelmann, 
& M. H. Soffel, Vol. 261, 183-188 



Masiero, J. R., Mainzer, A. K., Grav, T., Bauer, J. M., Wright, E. L., McMillan, R. S., 
Tholen, D. J., & Blain, A. W. 2012, The Astrophysical Journal, 749, 104 

Milani, A., Chesley, S. R., Sansaturio, M. E., Bernardi, F., Valsecchi, G. B., & Arratia, O. 
2009, Icarus, 203, 460 

Milani, A. & Gronchi, G. 2009, Theory of Orbit Determination (Cambridge University 
Press) 

Morbidelli, A. & Vokrouhlicky, D. 2003, Icarus, 163, 120 
Opeil, C, Consolmagno, G., & Britt, D. 2010, Icarus, 208, 449 

Ostro, S. J., Benner, L. A. M., Giorgini, J. D., Nolan, M. C, Hine, A. A., Howell, E. S., 
Margot, J. L., Magri, C, & Shepard, M. K. 2005, IAU Circular, 8627 

Pravec, P. & Harris, A. 2000, Icarus, 148, 12 

Pravec, P. & Harris, A. W. 2007, Icarus, 190, 250 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical 
Recipes in C (2nd ed.): The Art of Scientific Computing (New York, NY, USA: 
Cambridge University Press) 

Sekanina, Z. 1976, Icarus, 27, 265 

Shepard, M. K., Clark, B. E., Nolan, M. C, Benner, L. A. M., Ostro, S. J., Giorgini, J. D., 
Vilas, F., Jarvis, K., Lederer, S., Lim, L. F., McConnochie, T., Bell, J., Margot, 
J. L., Rivkin, A., Magri, C, Scheeres, D., & Pravec, P. 2008, Icarus, 193, 20 

Sitarski, G. 1992, The Astronomical Journal, 104, 1226 

— . 1998, Acta Astronomica, 48, 547 



-38 - 



Taylor, P. A., Margot, J. L., Vokrouhlicky, D., Scheeres, D. J., Pravec, P., Lowry, S. C, 
Fitzsimmons, A., Nolan, M. C, Ostro, S. J., Benner, L. A. M., Giorgini, J. D., & 
Magri, C. 2007, Science, 316, 274 

Vokrouhlicky, D., Chesley, S. R., & Matson, R. D. 2008, Astronomical Journal, 135, 2336 

Vokrouhlicky, D., Milani, A., & Chesley, S. R. 2000, Icarus, 148, 118 

Yeomans, D. K. 1991, The Astronomical Journal, 101, 1920 

Yeomans, D. K. 1992, The Astronomical Journal, 104, 1266 

Ziolkowski, K. 1983, in Asteroids, Comets, and Meteors, ed. C.-I. Lagerkvist, H. Rickman, 
171-174 



This manuscript was prepared with the A AS IATjrpC macros v5.2. 



-39 - 



■* 2 

X 





t- 


CO 


Tf 


LO 


■X' 




CO 


CO 


Th 


CN 


CN 


TP 


CI 


04 


TP 


o 


o 


CO 


o 


CO 


CN 


CN 


CO 


co 


CN 


co 


CO 


TP 


CO 


m 






TP 






CO 


TP 


CO 


CO 


CI 


CN 


Tf 


CM 


CM 


CI 


d 


d 


d 


d 


d 


d 


o 


o 


o 


o 


o 




d 


d 


d 


d 


d 


d 


d 


1 — 1 


TP 

I— I 


d 


d 


to 


00 




CO 


tp 


in 


CN 


in 


TP 


CO 


X 


CN 


o 


TP 


CM 






o 


o 


TP 


X 


in 


CN 


lO 


CN 




o 


04 




CO 


TP 


o 


X 


in 


TP 


CN 


04 


00 


m 


CO 


CO 


o 


CO 


CM 


o 


in 


o 


o 




TP 


04 


00 


lO 


CO 


CO 


I— I 


o 


04 


Ol 


00 


00 


X 


CO 


X 


00 




t— 


t- 


in 




TP 


CO 


CN 










































co 


o 


o 


t- 


lO 


I— I 


CO 




CM 


CO 


CI 


o 


CN 


TP 


TP 


X 


l> 


CN 


I-H 


I— I 


Ol 




00 


CN 


in 


o 




X 


I— I 


TP 


CM 


CD 


CM 


in 


04 


I-H 


o 


CO 


o 


CO 


TP 


i- 


TP 








I-H 


i-H 


CN 


CO 


CO 


TP 


1 — 1 


CM 


I— I 


I— I 


I-H 


I-H 




CN 


CO 


I— I 


CO 


I— I 


CN 


TP 


i-H 


TP 



1— ' "i. 



C3 



CM..CO.04 m .... t- 

d;;d:d:::;;;;;:t-i;::;cN 

co. .in. co tp . . . . Oi 

CM CM i-H in CM 



8 



l 

o 



E3 
< 



< 



a, 



CN 


o 


CO 




tp 




in 


o 


in 


o 


CO 


in 




■CO 


00 


■CO 


CO 


CN 


CO 


CM 


in 


CO 


CN 


d 


I-H 


d 


d 


CN 


d 


CO 


CM 


CN 


CN 


CO 


CM 


CM 


CN 


d 


TP 


CO 


CN 


TP 


1> 


i — i 


TP 


CO 


CO 


CO 


X 


tp 


in 


tp 


in 


TP 


t~ 


CN 


CN 


O 


i-H 


CO 




CM 


X 


X 


CN 


TP 


CO 


o 


X 


o 




rH 


H 




CN 




i— I 


in 


CO 


TP 


TP 


in 


CN 


co 


TP 


CO 


TP 


CO 


CI 


CO 


in 


Ol 






1 


1 


■ 


I 


— 

i 


1 




1 


I 




■ 








1 


1 


1 


i-H 
1 


■ 


■ 




CM 




CM 


o 


CN 


CN 


CM 


o 


C4 


Ol 




X 


X 


O 


X 


t- 


X 


CI 


CI 


CI 


CM 


Ol 


■CI 


















o 


o 


o 


o 


o 




o 


o 


o 


o 


o 


O 




o 


o 


o 


o 


O 


o 


O 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


CM 


CN 


CN 


CN 


CN 


CN 


CM 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CM 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CO 


CN 


N- 


■ 

X 


\h 




■ 

CO 


■ 

o 


■ 

CN 


1 

CI 


1 

l> 


00 




00 


CN 


X 


1 

lO 


1 

in 


in 


1 


TP 




d 


a: 


Ol 


in 


TP 


ID 


in 


I-- 


in 


X 


Tf 


t- 




t- 


t- 


X 


X 


t— 


t— 


in 


in 


CN 


Ol 


X 


X 


Ol 


Ol 


04 


Ol 


Ol 


Ol 


CI 


CI 


CI 


Ol 


Ol 


Ol 


Ol 


Ol 


Ol 


CI 


CI 


CI 


CI 


Ol 


Ol 


Ol 


lO 




lO 


1-H 


CO 


X 




X 




i-H 




i-H 


CN 


Ol 


X 






CO 






CO 


TP 


■CO 


CN 




CM 


CO 


CN 


in 




o 




m 




CO 


CN 


o 


i-H 






i-H 






CM 


i-H 


CM 


d 




d 


d 


d 


CN 

CO 




d 




d 




d 


d 


d 


d 






d 






CI 
CN 


d 


d 



o 


in 


m 


t- 


tP 


tP 


X 


CO 


o 


t- 


o 


X 


'CO 


o 




t~ 




'CO 






CO 


Ol 


t- 


o 


CM 


in 


t- 


CI 


Ol 




■CO' 


tP 


'CO 


CI 


CI 


■ m . 


CO 


O 


CM 

in 


o 

CO 


04 


d 

TP 


CM 
CM 
CM 


CO 
1 — 1 


CO 


CM 
CN 


TP 
i-H 


CO 

t- 


X 

d 


in 
t- 


o 
o 


t- 

CM 


CO 
CN 


1> 

d 


d 
i — i 


CN 
CM 


'CO 
CN 
TP 






d 




CO 












CO 


d 






i — i 









E3 
< 



X 


CN 

tP 
d 


o 
in 


o 

TP 


o 

i-H 


'CO 

m 


o 

CO 


o 

X 


O 

CO 


o 
q 


in 

CO 

i — i 


o 

CM 


o 

CM 


o 
q 


o 
in 


in 

'CO 

d 


4t 


o 

CO 


CN 

in 
d 


CO 
CI 

d 


CO 

q 


o 

CO 


o 
q 


d 


i-H 


CO 


i-H 


CN 


d 


in 


d 


i-H 


TP 


i-H 


i-H 


d 


q 

CN 


CN 


i-H 

CO 


d 


i-H 


CO 
X 


in 

CN 


in 

■CO 


X 
CO 


■CO' 

Ol 


tP 
CO 


in 

X 




Ol 
CO 


CO 
X 


CO 
tP 


'CO 
CO 


o 


'CO 

q 


X 


o 

m 


o 


CO 

t- 


t- 


Ol 
CO 


o 
1> 


CN 

TP 


X 


d 




d 


d 


00 


CO 
i-H 


in 


CM 

m 


CM 


CM 
CM 


d 


d 


d 


d 

CN 


CO 


CO 


Tf 

CO 


in 

i-H 


X 


TP 
i-H 


CO 
CN 


i-H 
i-H 


CN 


CN 
CM 


t- 
CN 


'CO 

in 


TP 
TT 


X 


tP 
■CO 


m 

TP 


Ol 
Tf 


CN 
CO 


CO 
X 


m 

CO 


m 


tP 


rH 

CN 


in 

TP 


o 
in 


o 

CO 


tP 

TP 


CM 
i — 1 


in 

CO 


CO 

in 


o 

TP 


'CO 

tP 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


CO 
TP 


Ol 


t- 

tP 


CO 


t- 

Ol 


in 

CM 


TP 
X 


o 

CM 


i-H 

CM 


X 

q 


X 

q 


X 
TP 


X 

q 


CM 
i-H 


CO 
00 


t- 

CM 


Ol 
CN 


CO 
X 


Tf 

q 


o 

CO 


CO 
CO 


CN 
t- 


■CO 

q 




d 


i-H 


i-H 


d 


i-H 


d 


CN 


i-H 


i-H 


i-H 


CN 


i-H 




i-H 


i-H 


i-H 


d 


i-H 


i-H 


CM 


i-H 





fa 



o 
ft 

<; 



o 
a. 



o 
J3 
— i 
cS 



^cj 

ft 







pq O g m < o 



O 



o 



J3 

PS 



Ph 



CO 

o 



CO 


■CO 


CN 


m 


CN 


o 


O 


o 


; — 1 


CO 


CO 


t- 


in 


Ol 


1~- 


t- 


CN 


o 


in 


CI 


CO 


CN 


CO 


'CO 


'CO 


X 


'CO 


CN 


tP 


X 


CO 


CO 


CO 


X 


'CO 


X 


m 


co 


o 


o 


TP 


CI 


CO 


o 




in 


X 


CO 


o 


'CO 


CO 


in 


CO 


in 


o 


X 


X 


Ol 


t- 


o 


— 


1 — I 


Ol 


CO 


o 


tP 




CN 


H 


r-l 


CN 


i-H 


CN 


i — i 


CO 


i — i 


CN 




i-H 


!-H 


co 


'CO 


CN 


CN 


t- 


t- 


i-H 


tP 




m 
























CO 










o 


CO 




iH 



i — I i-H (M 



-40- 



b- 


00 


CO 


X 






Oi 


Oi 


to 


CN 


b- 


CO 


o 


CN 


o 


00 


X 


TP 


1 — 1 


00 


to 




00 


(M 


to 


cs 


GO 


TP 


CO 


CO 


CO 


LO 


tP 


TP 


05 


q 


q 


q 




b- 


00 


q 


^* 


b- 


CO 


q 


d 


CN 


t-H 


T— I 


d 


CN 


I— I 


d 


i— i 


d 


d 


CN 








CN 


I— I 


CN 


i— i 


I— I 


CN 


T— I 




CC 


co 


to 


o 


m 


TP 


Oi 


to 


to 


CN 


o 


CN 


05 


00 


T— I 


00 


CO 


00 


b- 


to 


i— I 


'CO 


TP 


t? 


TP 


to 


q 


05 


00 


LO 


CN 


CN 


b- 


lO 


TP 


CO 


CO 


CN 


Oi 


Oi 


00 


00 


b- 


CO 


lO 


TP 


LO 


lo 


lo 


LO 


TP 


TP 


TP 


TP 


TP 


to 


to 


CO 


CO 


CO 


CO 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


Oi 


CN 


CO 


Oi 


co 


CN 


b- 


o 


LO 


TP 


LO 


00 


'CO 


TP 


LO 


05 


o 


CO 


TP 


to 




00 


co 


CN 


Oi 


Oi 




CN 




TP 


to 


to 


I— I 


CN 


CO 


TP 


o 


CN 


■00 


00 


to 


LO 


CN 


CN 


o 


■00' 



05 

d 



05 

TP 



'CO 

d 



l 

o 



b- 

d 



1= 
< 



I 

o 



c 


lO 


'CO 


CN 


05 


o 




Oi 




Oi 


o 


Oi 


b- 


CN 


LO 


■00' 


00 


o 


■* 




b- 


Oi 






o 


o 




O 






o 




o 




o 


o 




o 


o 


o 




o 




o 


o 




o 


o 


o 


O 


O 


o 


o 


o 


o 


o 


o 


o 


o 


O 


o 


o 


o 


o 


o 


c 


o 


o 


o 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 


CN 






05 


d 


LO 


CO 


■ 

Oi 


■ 

00 


■ 

CN 


■ 

CO 


■ 

CN 


CO 


CN 






CN 


■ 

lO 


4 


4 


■ 

LO 


■ 

LO 


CN 


CO 


X 


b- 


■00' 


CO 


lO 


b- 


00 


b- 


lO 


00 


00 


'CO 


00 


00 


1 o 


b- 


lO 


LO 


00 


lO 


lO 


LO 


b- 


Oi 


05 


05 


Oi 


Oi 


Ol 


Oi 


Ol 


Oi 


Oi 


Oi 


Oi 


Oi 


Oi 


Oi 


Ol 


Oi 


Oi 


Oi 


Oi 


Oi 


■Oi 


05 



a, 



to 
d 



CN 
TP 
CO 



I o 
CN 

TP 



O 

-a 



to 






o 






o 






CN 
















o 


O 








o 


to 


b- 


LO 


q 


CN 


Oi 


Tf 


LO 




■* 


LO 


'CO 


00 


CO 


CN 


b- 


CO 


00 


b- 


CO 


CO 


00 


TP 


d 


LO 


TP 


CN 


q 


■00' 




CO 




d 


CO 


TP 


■Oi 


LO 




q 


q 




d 


Tf 


^CH 




LO 




CN 


d 






CO 




d 






d 


CN 


d 


d 






CN 






CN 


CN 






to 


CO 


00 


b- 


o 


tH 


05 


lO 


o 


b- 


b- 


TP 


CN 


'CO 






CO 


b- 


CN 


CO 


CO 


CO 


LO 


TP 


q 


to 


b- 


CN 


■X' 


00 


q 


I— I 


i — i 


i— i 


q 


TP 


q 


00 




I— I 


q 


Oi 


b- 


LO 


CO 


TP 




00 


TP 


co 


CO 


oi 


00 


LO 


l> 


i — i 


b-^ 




TP 




d 


00 


LO 


CN 


d 


CN 


d 


LO 


d 




■to 












— 


CN 


i — i 




T— I 










CN 


i — I 




i— I 


CN 


■to 




CC 


'CO 


TP 


Tf 


b- 


T— I 


00 


b- 


TP 




00 


CO 


LO 


00 


b- 


rH 


TP 


b- 




TP 


lO 


'CO 


■to 


to 


b- 




I o 


LO 


LO 


TP 


i— i 


to 


^P 


TP 


LO 


CN 


o 


co 


LO 


to 


b- 




CO 


CN 


CN 


'CO 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


d 


05 


■Oi 


t- 


CO 


CN 


o 


CO 


b- 


o 


CO 


i— i 


TP 




1 o 


X 


Oi 


CO 


o 


Oi 


LO 


b- 


to 


CO 




b- 


CN 


CN 


CN 


o 


o 


to 


X 


o 


to 


'CO 


K 


■to 


00 


CN 


o 


to 


Oi 


CO 


CO 


CO 


1 o 








CN 






I— I 


i— I 


I— I 


i— i 


CN 


CN 








CN 


— 


CN 


o 


1 — 1 


i— I 




CN 



1 x 



a 
< 



O O £ 



o 

o 

m 



o 
z 

o 
o 
o 



fa 

a 

CN 
O 

o 



s 



CN 
O 

o 



to 





f3 



o 


o 


CN 


H 


CO 


■to 


Oi 


b- 


CO 


tP 


CO 


05 


LO 




Oi 


CN 


O 


b- 


CN 


LO 


O 


o 


Oi 


CC 


CO 


o 


CO 


1 o 


LO 


CO 


Tf 


Oi 


to 


to 


LO 


TP 




o 


o 


LO 


05 


CD 


to 


LO 


■Oi 


b- 


CC 


'CO 


CO 


o 


b- 


b- 


b- 


Oi 


LO 


O 


to 


o 




Oi 




CN 


Oi 


i— I 


to 


Oi 


to 


'CO 






1 o 


o 


CN 


CN 


CO 


tP 


TP 


I— I 


TP 


b- 


CN 


LO 


00 


00 


CN 


X 




to 


Oi 


X 


TP 












CN 








TP 
CN 










■to 






CO 






Oi 


CO 


LO 





41 - 



2 



CO 




I 

o 



o 
O 



co 
o 

3 



a. 



3 



3 



< 

E 



o oi t- 

r-i 00 



to 
d 



Ol Ol Ol Ol Ol oi 



CN 








TP 




o 




t- 








TP 


Ol 


■CO' 




lO 








CO 


00 


lO 




00 








CO 


CO 


ci 












CO 








o 












00 






o 


00 


00 


Ir- 


TP 


TP 




CO 




CO 




o 




CN 






d 


d 






CN 






o 


CO 


CN 


CO 


ID 


Ol 


CO 


CN 


CN 






TP 


cn 


CN 


■CO' 


1C 


c4 




CO 


CO 




lO 


CO 




(M 






CO 






CN 








00 


t~- 




Oi 


GO 


CO 


CO 


TP 


TP 


CN 


CO 


TP 


CN 


CD 


o 


o 


o 


o 


o 


o 


o 


o 


CO 


00 




H 


TP 


CO 


t- 


00 




TP 


TP 


o 


'CO 1 


Ol 


Ol 





X 

> 



a 



fa 





o 


CO 


co 


t-H 


CO 


TP 


TP 


CO 


CO 




CO 


lO 


CN 


lO 


lO 


'X 


CO 






CO 


Ol 


lO 


CN 




TP 




o 


■CO 


CN 


■CO 


00 








rH 


T— I 






00 




I— I 































-42 - 



Table 4: Statistical properties of observed Yarkovsky rates and efficiencies. 









&bs(da/dt) x 10~ 4 AU/Myr 




Yarkovsky rate 




mean 


weighted mean 


median 


stdev 


objects with f Y <2x 10" 


-5 


7.6 


4.4 


5.6 


6.4 


objects with f Y >2x 10" 


-5 


27.0 


18.5 


20.1 


18.7 


all objects 




10.4 


5.2 


7.3 


11.4 








fy X 10- 


-5 




Yarkovsky efficiency 




mean 


weighted mean 


median 


stdev 


objects with f Y < 2 x 10" 


-5 


0.67 


0.53 


0.50 


0.51 


objects with f Y > 2 x 10" 


-5 


4.50 


7.47 


3.01 


3.38 


all objects 




1.22 


0.89 


0.65 


1.91 



